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ABSTRACT 



Context. Size measurements of young star clusters are valuable tools to put constraints on the formation and early dynamical evolution 
of star clusters. 

Aims. We use HST/ACS observations of the spiral galaxy M51 in F435W, F555W and F814W to select a large sample of star clusters 
with accurate effective radius measurements in an area covering the complete disc of M5 1 . We present the dataset and study the radius 
distribution and relations between radius, colour, arm/interarm region, galactocentric distance, mass and age. 

Methods. We select a sample of 7698 (F435W), 6846 (F555W) and 5024 (F814W) slightly resolved clusters and derive their effective 
radii (R e g) by fitting the spatial profiles with analytical models convolved with the point spread function. The radii of 1284 clusters 
are studied in detail. 

Results. We find cluster radii between 0.5 and ~10 pc, and one exceptionally large cluster candidate with R e g = 21.6 pc. The median 
i? e ff is 2. 1 pc. We find 70 clusters in our sample which have colours consistent with being old GC candidates and we find 6 new "faint 
fuzzy" clusters in, or projected onto, the disc of M51. The radius distribution can not be fitted with a power law similar to the one for 
star-forming clouds. We find an increase in R e g with colour as well as a higher fraction of clusters with B—V > 0.05 in the interarm 
regions. We find a correlation between /? c(f and galactocentric distance (R G ) of the form R cS oc i^ 12±0 - 02 , which is considerably weaker 
than the observed correlation for old Milky Way GCs. We find weak relations between cluster luminosity and radius: R^ oc £ 015±002 
for the interarm regions and R eff oc Zr - 1 1±0 01 for the spiral arm regions, but we do not observe a correlation between cluster mass and 
radius. 

Conclusions. The observed radius distribution indicates that shortly after the formation of the clusters from a fractal gas, the radii of 
the clusters have changed in a non-uniform way. We find tentative evidence suggesting that clusters in spiral arms are more compact. 

Key words, galaxies: individual: M51 - galaxies: star clusters 



1. Introduction 

One of the most striking questions related to star cluster forma- 
tion concerns the transition from the densest parts of the star- 
forming Giant Molecular Clouds (GMCs) to the clusters that 
emerge from them. Observations show that the masses and radii 

1 /2 

of cl ouds follow a clear relation of the form Rcmc k ^gmc 
(e.g. lLarsonlll98lt ISolomon et al.lll987l) down to the scale of 
the star forming clumps (Willia ms et al.l [l995l) . Such a rela- 
tion is expected for self-gravitating clouds, which are in virial 
equilibrium and pressure bounded by their surroundings (e.g. 
Elmegree^[l98§. The same mass-radius relation is observed 
for complexes of multiple clusters dElmegreen & SalzeJ 1999; 
Elmegree n et alJl200U iBastian et al.ll2005ab . Elliptical galaxies 
and very massive (> 1 6 5 M B ) stellar clus ters follow a simi- 
lar relation (R oc M 62 (Hasegan et al. 2005)), which advocates 
a possibility of f orming very massive clusters by merging of 
low m ass clusters (Fellh auer & Kr oupa 2005; Kissler-Pati g et al.l 
2006). 



However, for individual star clusters that emerge from the 
star forming cl ouds/clumps, a relation between mass and radius 
is not presen t dVanden Berg h et all [1991 IBastian et alJl20 05b: 
I jordan et"aL 2005b or at le ast strongly weakened dZepf et al 
Il999t iHunter et al.l 120031; [Mackey & Gilmord l2003t Larser 



20041 iLee etai1l2005l) . Since this is already the case for very 
young clusters, it indicates that during, or shortly after (< 
10 Myr) the transition from clouds to clusters the mass and/or 
the radius of the objects change. 

These changes in mass and/or radius are likely to be 
reflected in changes in the mass and radius distributions 
d Ashman & ZepflfeOOll) . On the one hand, however, the mass 
distributions of both clouds and clusters show great similar- 
ities. Both can be approximated by power laws of the form 
N(M) dM oc M~ a dM, with the index a in the range of 1.5-2.0 
(see [Senders et alJI 19851; ISolomon et all 19871; lHarris & Pudritzl 
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1 994 :lBrand & Wouterlootll 19951; lElmegreen & Falgarondl99l 
Fukui et a ll 1200 ll for c l ouds and see e.g . [Zhang & Fall 1999; 
Ashman & Zepfl hooU iBik et all [20031 Ide Griis et al] [2003; 
Hunter et alj 120031 for clusters). Recently, some studies have 
found evide nce for an upper ma s s truncation of the clu ster mass 
distribution iGieles et all d2006al) ; iDowell et all d2006l), which is 
also found for the mass distributions of GMCs (Rosol owskvl 
2005). On the other hand, the radius distributions are less well 
constrained, especially for extra-galactic star clusters. If we 
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approximate the radius distributions by a power law of the 
form N(R)dR oc R n dR, the average value of the index for 
GMCs is observed to be n_=_ 3.3 ± 0.3 dHarris & Pudritzlll994l: 
Elmegreen & Falga rond[l996h . consistent with the gas having 
a fractal structure with a fractal dimension n - 1 = 2.3 + 0.3 
( Mandelbrot! 1 1 983l lElmegreen & Elmegreen! 2001b . For young 
clusters in NGC 3256 IZepf et all d!999h find 77 « 3.4, while 
Bastian et al. (2005, from now on referred to as "B05") find 
77 = 2.2 ± 0.2 for star clusters in the disc of M51. This differ- 
ence, however, seems to be ca used by the erroneous a ddition of 
+ 1 to the index in the result of lAshman & ZepJ^OOlb Pl 

Our understanding of whether or not the mass and radius 
distribution of clouds and clusters are similar, and to which de- 
gree, is directly coupled to our understanding of star formation 
and the early evolution of star clusters. Besides that, possible 
explanations for the lack of a mass-radius relation for clusters 
which change the mass and/or the radius of the clusters in a non- 
uniform way, are likely to affect the mass and radius distributions 
(Ashman & Zepf 2 0011) . It is therefore important to put better 
constraints on these distributions, and in this work we will focus 
on the radius distribution of young star clusters. 

In the study presented here we exploit the superb resolution 
and large field-of-view of the new HST Advanced Camera for 
Surveys (ACS) observations of M51, taken as part of the Hubble 
Heritage Project. These observations allow us to measure the 
radii of a large sample of clusters in an area covering, for the first 
time, the complete disc of M51 and its companion, NGC 5195 
at a 0.05" resolution. In this work, which is the first in a series 
of papers, we present the dataset and we study the radius distri- 
bution for the complete cluster sample and for cluster samples 
with a different background surface brightness ("background re- 
gions"). The differences in background regions are likely to re- 
flect differences in environmental conditions, which could have 
an impact on the early evolution of star clusters. 

The rad ii of star clusters in M5 1 have already been studied 
by B05 and lLee et alj 120 05). However, these studies used lower 
resolution HST/WFPC2 data and were not covering the complete 
disc. Besides this, we use different selection criteria for the clus- 
ters than B05, based on the clusters actually being resolved and 
clearly separated from nearby contaminating sources. In combi- 
nation with the larger field-of-view and the higher resolution of 
the ACS data, this leads to a larger sample of clusters, divided in 
different background regions. 

The structure of this paper is as follows: in §[2] we describe 
the dataset, the selection of sources and background regions and 
the photometry. The radius measurements are described in § [3] 
and in §|4]we describe experiments with artificial clusters to de- 
termine the accuracy and detection limits of our methods. Our 
selection criteria and a description of our cluster sample are 
presented in § [5] followed by a comparison between ACS and 
WFPC2 data in § |6] The radius distribution is presented in § |7j 
and we search for correlations between mass, radius and galac- 
tocentric distance in §[8] In ^9] we finish with the summary and 
conclusions. 



Table 1. Overview of the filters and exposure times used. The 
limiting magnitude is for a point source with S/N = 5. 



Filter 




Exposure 


time 


Limiting magnitude 


F435W 


4 


x 680 s = 


2720 s 


27.3 m B 


F555W 


4 


x 340 s = 


1360 s 


26.5 my 


F814W 


4 


x 340 s = 


1360 s 


25.8 ni! 


F658N+N[II] 


4 


x 680 s = 


2720 s 





1 Fitting a function of the form \og(N(r)) = a + b ■ log(r) results in 
b = -7]. However, using logarithmic binning, one fits log(A'(log(r))) = 
a + b ■ log(r), in which b = 1 - 77. This extra term +1 can easily cause 
confusion when comparing di fferent distributions. Also see B05 and 
Elmegreen & Falgarone ( 1996). 



2. Observations, source selection and photometry 

2.1. Observations 

We make use of the HST I ACS dataset of M5 1 (NGC 5 194, a late 
type Sbc galaxy), taken as part of the Hubble Heritage Project in 
January 2005 (proposal ID 10452, PI: S. V. W. Beckwith). The 
dataset consists of 6 ACS pointings using the Wide Field Channel 
(WFQ in F435W (~ B), F555W (~ V), F814W (~ /) and F658N 
(~ Ha), with 4 dithered exposures per filter, as is summarized 
in Table Q] The obs ervations were reduc ed and drizzled into one 
mosaic image bv lMutchleretalJ d2005l) . In summary, the stan- 
dard ACS pipeline (CALACS) was used for bias, dark and flat- 
field corrections of the individual dithered images. The corrected 
images were then com bined into one mosaic image using the 
Mult iDrizzle task dFruchter & Hookll2002l: iKoekemoer et alj 
|2002|) . which also corrects for filter-specific geometric distor- 
tion, cosmic rays and b ad pixels. For a comp lete description of 
the dataset we refer to Mut chler et al.l d2005l) and the M5 1 mo- 
saic website (http://archive.stsci.edu/prepds/m51/). 
For details on the sta ndard pipeline process ing we refer to the 
ACS Data Handbook (Pavlovsky et al. 2005). 

The resulting mosaic image covers a region of 430 x 610" 
(= 17.5 x 24.8 kpc) with a resolution of 2.0 pc per pixel, 
where we assumed a d istance modulus of m - M = 29.62 from 
iFeldmeier et alJd 19971) . i.e. a distance of 8.4 + 0.6 Mpc. The cov- 
ered region is large enough to include the complete spiral disc of 
M51, as well as its companion NGC 5195 (a dwarf barred spiral 
of early type SB0), while at the same time the resolution is good 
enough to resolve stellar cluster candidates, i.e. to distinguish 
them from stars by measuring their sizes and comparing these to 
the Point Spread Function (PSF) of the HST '/ACS camera. 

2.2. Source selection 

For source selection we used the SExtractor package 
dBertin & Arnoutslll996l version 2.3.2). SExtractor first gener- 
ates a background map by computing the mean and standard de- 
viation of every section of the image with a user defined grid 
size for which we choose 10 x 10 pixels. In every section the 
local background histogram is clipped iteratively until every re- 
maining pixel value is within +3cr of the median value. The 
mean of the clipped histogram is then taken as the local back- 
ground value. Every area of at least three adjacent pixels that 
exceeded the background by at least 5cr was called a source. For 
details on the background estimation and the source selection 
w e refer to the SExt ractor user manual (Be rtin & Arnoutsl fl996) 
or iHolwerdal (|2005). The coordinates of the sources in F435W 
F555W and F814W were matched and only sources that were 
detected in all three filters within two pixel uncertainty were 
kept. This resulted in a list of 75 436 sources, including clus- 
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ter candidates but also many stars and background galaxies. We 
did not apply any selection criteria based on shape, sharpness or 
size during the source selection with SExtractor. However, in §|5] 
we use individual radii measurements to select a large sample of 
cluster candidates from the source list. 

2.3. Background selection 

In order to study any possible relations between radius and en- 
vironment, and to perform completeness and accuracy tests for 
different background levels, we divided the image in three re- 
gions according to the background surface brightness. These 
background regions were selected by smoothing the F555W im- 
age with a boxcar average of 200 pixels. Consequently, pixels 
with a value < 0.19 e~/s (corresponding to a surface bright- 
ness fi > 21.02 mag arcsec -2 ) were flagged "low background". 
Pixel values > 0.19 e"/s and < 0.25 e"/s (20.72 < fi < 
21.02 mag arcsec -2 ) were flagged "intermediate background" 
and pixels with a value > 0.25 e~/s (jj. < 20.72 mag arcsec -2 ) 
were flagged "high background". These values were chosen be- 
cause they resulted in a contour map, shown in Fig.Q] in which 
the high background region clearly follows most of the inner 
spiral arms, covering most areas that would be marked "high 
background" if the selection would take place by eye. The "in- 
termediate" region should be considered as a transition region 
to clearly separate cluster samples within a low and high back- 
ground region. 

2.4. Point spread function 

For our photometry, radius measurements and artificial cluster 
experiments we need a filter dependent PSF. Since there are not 
enough bright, isolated point sources in the M5 1 mosaic image 
to determine the PSF accurately, the PSF was observationally 
determined from a crowded star field on a drizzled image of 
the Galactic globular cluster 47 Tuc (NGC 104, HST proposal 
GO-9281, PI: J. E. Grindlay). For each filter a separate PSF was 
determined. 

With drizzled data, the size of the PSF depends on the choice 
of the drizzle kernel and the accuracy with which the flux from 
multiple images is combined. We note that the image of 47 Tuc 
was drizzled in a slightly different way than the M51 image. 
The 47 Tuc images were drizzled usi ng a square kernel with a 
size of one pixel (i.e. pixfrac = 1.0 dFruchter & Hookll2002lV ). 
while a Gaussian kernel with pixfrac = 0.9 was used for M5 1 . 
Therefore, we expect some differences between the PSFs, with 
the M51 PSF possibly being larger than the 47 Tuc PSF. This 
could le ad to an overe stimation of the measured radii. However, 
tests by Larseri] d2004l) have shown that the lower limit down to 
which Ishape can detect a source as being resolved is ~10% 
of the FWHM of the PSF. At the distance of M51 and with a 
FWHM of the ACS PSF of -2.1 pixels, this corresponds to an 
effective radius (/? e ff) of ~0.5 pc. We indeed find a very strong 
peak in the R e g distribution of all the measured sources below 
0.5 pc, consistent with the majority of the sources (faint stars) 
being fit as point sources. 

This peak of point sources implies that the measured radii are 
not biased towards larger values. We therefore conclude that the 
empirical PSF we use, although drizzled in a slightly different 
way than the image of M51, is not too small. It shows that other 
effects on the PSF, like combining the flux of multiple separate 
images, are more important than the differences in the drizzle 
parameters. In §|6]we also show that there are no large systematic 




Fig. 1. The contours outlining the three background regions, su- 
perimposed on the HST '/ACS mosaic image in the F555W pass- 
band. The solid lines enclose "high background" regions, and ev- 
erything outside the dashed lines indicates the "low background" 
regions. The regions in between solid and dashed lines are called 
"intermediate background". 



offsets between our measured radii and the radii of some clusters 
that were independently measured by B05 using WFPC2 data. 
We can therefore use the empirical PSF from the image of 47 Tuc 
and we will use as a lower-limit for the measured /? e ff a value of 
0.5 pc. 

2.5. Photometry 

We performed aperture photometry on all the sources in the 
source list using the IRAF^DAOPHOT package. We used a 5 
pixel aperture radius and a background annulus with an inner 
radius of 10 pixels and a width of 3 pixels. 

The aperture correction (AC ()5 ») for resolved sources from 
the 5 pixel aperture to 10 pixels (= 0.5") depends on the size 
of the source. Larger sources will have more flux outside the 
measuring aperture, and therefore need a larger (i.e. more nega- 
tive) aperture correction. We measured the aperture corrections 



2 The Image Reduction and Analysis Facility (IRAF) is distributed by 
the National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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on artificial sources with different effective (i.e. projected half- 
light) radi i between ~ 1 pc and ~ 5 pc, generated by the BAOlab 
package ( ILarsenl 1 19991 120041). For these artificial sources we 
used Moffat profiles ((Moffat 1969]) with a power-law index of 
-1.5, which we convolved with the filter dependent PSF of the 
HST/ACS. The aperture correction was then measured by com- 
paring the photometry between a 5 and 10 pixel aperture. 

The measured aperture corrections in all the three filters 
(F435W, F555W and F814W) show a clear linear relation with 
the measured size of the analytical cluster. Fig. [2] shows the re- 



F435W 



suits for F435W. The relations between AC 



0.5 



and measured 



effective radius (Reg) for the different filters can be written as: 



F435W : 



F555W : 



AC, 



0.5" 



mag 
AC n 



-0.045 - 0.043-1 — 
pc 



"0.5 



F814W : 



mag 
AC n 



-0.045 - 0.042 • 



eff 



-0.5 



mag 



pc 

-0.057 - 0.043 ■ | — 
pc 



(1) 



For a source which has a typical measured size of 3 pc (see §[7]), 
this corresponds to an aperture correction of -0.17, -0.17 and 
-0.19 mag in F435W, F555W and F814W, respectively. 

We could use Eq. Q]to apply a size-dependent aperture cor- 
rection to every source based on their radius measurement. 
However, this could introduce new and unknown systematic un- 
certainties due to the limitations of the radius measurements. We 
therefore decided to use a fixed aperture correction, determined 
for a 3 pc source. We are aware that this introduces uncertainties 
in the flux as a function of the sizes of the sources ( Ander set al.l 
2006). We will underestimate the flux with ~0.3 mag for a 9 pc 
source and overestimate the flux with ~0.1 mag for a 0.5 pc 
source (also see §[4]i. However, it is preferred to deal with these 
systematic uncertainties of known origin rather than introducing 
new uncertainties due to the less well determined uncertainties 
of the sizes. Moreover, since the coefficients in Eq. Q]are very 
similar for the different filters, uncertainties in the aperture cor- 
rections practically cancel out when we consider colours instead 
of fluxes. Nevertheless, when necessary we will mention how 
our results change if we use the size-dependent aperture correc- 
tions. 

The filter dependent aperture corrections from 0.5" to in- 
finity (the "infinity corrections") were taken from Table 5 of 
Sirianni et al. (2005), which were determined for point sources. 
In principle, with these infinity corrections for point sources we 
are slightly underestimating the infinity corrections for resolved 
sources. However, in §|4]we will show, using artificial cluster 
experiments, that with the infinity corrections for point sources 
we are not introducing systematic offsets in the photometry for 
3 pc sources. 

A final correction for Galactic foreground extinction of 
E(B-V) = 0.038 in the direction of M5 1 was applied, according 
to Appendix B of Schle geTet al.l l[l998). This corresponds to an 
additional correction in F435W, F555W and F814W of -0.152, 
-0.1 17 and -0.068 mag, respectively. 

We did not apply any CTE corrections, since on the drizzled 
mosaic image the exact location of every source on the CCD is 
not easily retrieved, nor did we do photometry on the single (i.e. 
"un-stacked") exposures, both of which are necessary to calcu- 
late the CTE corrections. We have estimated the CTE corrections 
to be of the order of -0.02 mag and therefore ignoring them does 
not lead to large systematic effects. We also did not apply a cor- 
rection for the impact of the red halo (Si rianni et al.l l200~5b on 




Fig. 2. The aperture correction (AC„ -//) for the flux between a 5 
and 10 pixel (= 0.5") aperture in F435W, versus the effective ra- 
dius R e ff of analytical sources at the distance of M5 1 (diamonds). 
The dashed line is a linear fit to the data of the form of Eq. Q] 



our F814W photometry, since clusters in the spiral disc of M5 1 
are mainly blue objects and the red halo effect is most significant 
for very red o bjects observed in th e F850LP filter. Using Tables 
6 and 8 from lSirianni et all (2005) we estimate that the error in 
the infinity correction for our F814W photometry would be of 
the order of ~ 0.01 mag, if the clusters would be red objects. 
This shows that the red halo effect has no significant effect on 
our photometry. 

3. Radius measurements 

We exploit the resolution of the ACS camera by measuring the 
effective radii of all the 75 436 detected sources in F435W, 
F555W and F814W. These radii can than be used to distinguish 
slightly resolved stellar cluster candidates from stars (§ [5) and 
to study the size distribution of a large sample of stellar cluster 
candidates (§ |7j». With "slightly resolved" we mean that the clus- 
ter candidates have an intrinsic size which is comparable to or 
smaller than the FWHM of the PSF. 

For the radius measurements we use d the I shape r outine , 
which is part of the BAOlab package (ILarsenl 1 19991 |2004|) . 
I shape convolves analytic profiles for the surface brightness 
distribution of a cluster with different effective radii with the 
PSF and then fits these to each source in the data. The best fitting 
R e g is then determined by minimizing the x 1 m an iterative pro- 
cess. For the analytic profiles we used the same ones as for the 
measured aperture corrections described in § 12.51 namely Moffat 
profiles with a power-law index of -1 .5 (i.e. a Moffat 15 profile). 
These profiles were found to be the best-fittin g profiles to young 
stellar clusters in the LMC (lElson et al.lll987l) . 

Because the M51 ACS data was drizzled, the cores of the 
surface brightness profiles of the young clusters in M51 could 
have been slightly changed. We did not quantify this effect, but 
instead stick to the Moffat 15 profiles, since the outer regions of 
the profiles, which in the case of Moffat 15 profiles approximate 
power laws, are not expected to change. 

The average surface brightness p rofile of Galactic globula r 
clusters (GCs) i s a Ki ng 30 profile dKindll962t |Harrislll996T) . 
Tests by ILarsenl (1 19991) have shown that when clusters that fol- 
low a King profile are measured using a Moffat 15 profile, the 
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effective radius is reproduced quite well. Therefore, even in the 
case where the true profiles of stellar clusters in M5 1 are not per- 
fect Moffat 15 profiles, the use of this profile will likely not lead 
to large systematic errors. 

The radii of the sources were measured using the flux within 
a 5 pixel radius around the centre of the source (i.e. using an 
I shape fitting radius of 5 pixels). To avoid neighbouring sources 
to affect the radius measurements, we rejected all sources which 
have a neighbour in the original source list within 5 pixels (see 
§[5]for a full description of the applied selection criteria). 

Since stellar cluster profiles are almost never perfectly spher- 
ically symmetric, we fitted them with elliptical profiles. We 
transformed the measured FWHM (in pixels) along the major 
axis into an effective radius according to the formulae from the 
I shape user's guide: 

R eB ~ 0.5(1 + ratio) x 1.13 x 2.036 • FWHM, (2) 

which gives R e ff in parsecs. The factor 2.036 accounts for the 
size of a pixel in parsecs at the distance of M5 1 , the factor 1 . 1 3 is 
the conversion factor from FWHM to R e s for a Moffat 15 profile 
and the term 0.5(1 + ratio) is a correction for the elliptical profile 
(using the derived aspect ratio). Because of the correction for 
ellipticity, Eq.|2]gives the radius of a circular aperture containing 
half the total light of an elliptical profile. This way we have a 
single value for the effective radius of elliptical sources and we 
still preserve information about the aspect ratios and position 
angles of the sources for future studies. 

4. Artificial cluster experiments 

To test how our selection of stellar cluster candidates (§[5]l de- 
pends on the brightness and size of the cluster and the back- 
ground region, and how accurate our radius measurements are, 
we performed tests with artificial clusters for all three filters, dif- 
ferent background regions and different cluster sizes. The results 
of these tests will be used in §|5]to select a sample of stellar clus- 
ter candidates with accurate radii. 

First we created artificial clusters using the MKCMPPSF 
and MKSYNT H tasks, which are part of the BAOlab package 
dLarsen|["l999l) . For every filter we convolved the PSF with a 
Moffat profile with a power-law index of - 1 .5 and effective radii 
between 1 and 9 pc, in steps of 2 pc. These artificial clusters 
were then scaled to a range of magnitudes between 18 and 26 
mag with steps of 0.5 mag. For every magnitude 100 clusters 
were added at random locations to every background region on 
the mosaic image by combining MKSYNTH with the imarith 
task in IRAF. We made sure that the minimum distance between 
these random locations was at least 30 pixels, in order not to 
introduce artificial crowding effects. 

We then performed the cluster selection on these sources 
in a similar way as with the normal data. We performed these 
tests for every filter individually, since taking into account the 
matching of every source in three filters, as we did with the nor- 
mal data, would imply creating artificial clusters with a range 
of colours (i.e. ages) for every magnitude, drastically increasing 
the computing time. However, by comparing analytical spectral 
energy distributions (SEDs) from GAL EV simple stellar popula- 
tion (SSP) models ( Schu lzet al J2002I) to the measured detection 
limits, we tested that for cluster ages up to ~ 1 Gyr, the F814W 
filter is always the most limiting filter for detecting clusters. We 
therefore do not expect that using the results of these tests for in- 
dividual filters is introducing large age-dependent biases in the 
derived detection limits. 




18 20 22 24 26 

Input F435W Mag. 



Fig. 3. The recovered fraction of artificial clusters as a function 
of magnitude and R e ff for the three background regions on the 
F435W image (i.e. "completeness profiles"). The horizontal dot- 
ted line indicates where 90% completeness is reached. The mag- 
nitudes at which the completeness profiles reach this 90% com- 
pleteness limit are summarized in Fig. [4] 



The artificial clusters were recovered by running SExtractor, 
followed by photometry in all three filters and by perform- 
ing size measurements in F435W and F555W. For F435W and 
F555W, we considered a cluster to be recovered when we found 
a resolved source (i.e. with a R e ff larger than our adopted lower 
limit of 0.5 pc (see § (5j and a y 1 which is lower than the x 1 
using only a PSF) within 1 pixel from the input coordinate and 
with a distance to the nearest neighbouring source of >5 pixels. 
For F814W we considered a cluster to be recovered when we 
found a source within 1 pixel from the input coordinate and with 
a distance to the nearest neighbouring source of >5 pixels. The 
recovered fraction as a function of F435W magnitude for the dif- 
ferent background regions and input radii is plotted in Fig. [3] For 
the F555W and F814W passbands the results are similar, except 
that the turn-off of the completeness curves happens at brighter 
magnitudes due to a lower S/N ratio of the F555W and F814W 
data (see below). 

The recovered fraction shown in Fig.[3]is scaled to the num- 
ber of clusters recovered at magnitude 18. We note, however, 
that a certain fraction of even the brightest artificial clusters is 
not recovered due to their vicinity within 5 pixels from a neigh- 
bouring source. Our completeness tests show that these initial 
losses will be ~3, ~13 and ~28% for the low, intermediate and 
high background, respectively. This shows that due to crowding 
effects one can never select a sample which is 100% complete. 
Depending on the selection criteria and background region, one 
will lose up to ~28% of the initial sample present in the data. 
This number will be even higher when one considers that young 
clusters are not randomly distributed across the spiral disc, but 
will mostly be clustered themselves in high background regions. 

Fig. [3] shows that larger clusters are harder to recover than 
smaller clusters with the same brightness. This is expected, since 
larger clusters have a lower surface brightness, which makes 
them easier to blend into the background. We define the mag- 
nitude at which 90% of the artificial cluster was recovered to 
be the 90% completeness limit. The values we found this way 
for the different filters, cluster sizes and background regions are 
plotted in Fig. El 
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Fig. 4. The 90% completeness limits versus effective radius 
(i? eff ), determined for F435W, F555W and F814W for the three 
background regions. 

For our current study, where we look at the radius distribu- 
tion of stellar clusters, it is not only important to detect clusters 
by measuring their radii, but the measured radii also have to be 
accurate. To test the limitations of our radius measurements, we 
looked at how the difference between input and measured radius 
of the artificial clusters depends on magnitude, input size and 
background region. In Fig. [5] we plot the 50th (i.e. the median), 
68th and 90th percentile of 6 versus the magnitude in F435W, 
where we define 5 as the relative difference between input and 
measured radius: 



6 = 



l^eff.in - ^eff,( 



R 



(3) 



eff,in 



We note that the pth percentile is the value such that p% of the 
observations (5) is less than this value. Fig. [5] shows for exam- 
ple, that when we select cluster candidates with F435W < 22.6 
mag, ~68% of the sample will have a radius uncertainty smaller 
than 20%. The figure also indicates that the radius measure- 
ments are the most limiting factor in the detection of stellar clus- 
ter candidates: at the 90% completeness limit for a 3 pc source 
(F435W = 24.2 mag), about 50% of the recovered sources are 
likely to have inaccuracies in their radii larger than 25^0%. 
Therefore, in §[5] we will select magnitude limits in F435W and 
F555W brighter than the 90% completeness limits of these pass- 
bands. Since we will not use the radius measurements in the 
F814W passband, in this passband the 90% completeness limit 
will be used in the selection of the sample. 

We also used the results of our artificial cluster experiments 
to test the robustness of our photometry and the accuracy of the 
applied aperture corrections described in § 12.51 In Fig.|6]we show 
the mean difference between the measured magnitude and the 
input magnitude (A mag) versus F435W magnitude for differ- 
ent sizes in the high background region. We applied the constant 
aperture correction for a 3 pc source accor ding to Eq.[T] a nd the 
infinity correction for point sources from Siria nni et alj {2005) 
that we described in § 12.51 First of all, Fig. [6] shows that the ap- 
plied aperture correction is very accurate, since the photometry 
of 3 pc clusters is almost perfectly reproduced to ~22.5 mag. 
This shows that applying the point source infinity correction to 
3 pc sources does not introduce systematic offsets in the photom- 
etry. Second, the range covered in A mag shows that when there 



Fig. 5. The median (top), 68th percentile (middle) and 90th per- 
centile (bottom) of 6 (see Eq.[3]l versus F435W magnitude for the 
artificial clusters with different radii in the high background re- 
gion. Below the horizontal line the uncertainty of the radii mea- 
surements is smaller than 20%, which is achieved for more than 
68% of the clusters up to 9 pc brighter than F435W = 22.6 mag 
(indicated by the arrow). 
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Fig. 6. The mean difference between the measured magnitude 
and the input magnitude of the recovered artificial clusters ver- 
sus F435W magnitude for different sizes in the high background 
region. The truncated line for the largest and faintest clusters 
means none of these artificial clusters were recovered. 

is no information about the radius of the cluster, the uncertainty 
in the photometry can be as large as ~ 0.5 mag for clusters with 
radii in the range 1-9 pc. 



5. Selection of the sample 

We used the radius estimates to distinguish the resolved clusters 
from unresolved objects. In this section we select two different 
cluster samples: a "resolved sample" with clearly resolved clus- 
ters, and a smaller subset from this sample, a "radius sample". 
The radius sample satisfies extra criteria that make the radii more 
reliable, and will be used to study the radius distribution (§ [7]i 
and the correlations between mass, radius and distance (§ [8). 
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The resolved sample will be used in anoth er study of the lumi- 
nosity function of stellar clusters in M51 (Gie les et al.l l2006a; 
lHaas et al.|[2007l in prep.). 

The cluster selection process is hampered by various fac- 
tors like an irregular background (spiral arms, dust lanes), con- 
taminating background galaxies and crowding effects which 
causes many sources to be (partially) blended with neighbouring 
sources, biasing the radius measurements. We tried to automate 
the selection of stellar cluster candidates as much as possible, 
taking into account all these factors. However, it was unavoid- 
able to subject the automatically selected sample to visual in- 
spection, to filter out any remaining contaminants or the stellar 
cluster candidates of which the radii measurements could not be 
trusted. 

5. 1 . Selection criteria 

1. Our first two selection criteria were concerned with the 
sources actually being resolved. As mentioned in §[3] we will 
use as a lower-limit for the measured R e g a value of 0.5 pc, 
since below this radius we observe a strong peak of unre- 
solved sources (most likely bright stars). Therefore, our first 
criterion for the selection of stellar cluster candidates was to 
select sources with /? e ft > 0.5 pc. We applied this criterion to 
the radius measurements in both F435W and F555W. We did 
not apply a criterion to the measured radii in F814W, since 
that would have restricted our sample too much due to the 
lower signal-to-noise ratio in this passband. 

2. Not only should a stellar cluster candidate be resolved when 
we fit a Moffat 1 5 profile, but using the profile convolved with 
the PSF should also result in a better fit than fitting the can- 
didate with the PSF alone. A second criterion was therefore 
to use the x 1 of the radius fit using the Moffatl5 profile, 
which should be lower than the ^ of a fit using only the 
PSF:^/4<1. 

After these criteria there were still many contaminants in our 
sample, e.g. patches of high background regions in between dark 
dust lanes that were detected as a source, background galaxies, 
blended sources and crowded regions. The next criteria were 
used to remove contaminants and to select only cluster candi- 
dates with accurate radii: 

3. Following the results of our artificial cluster experiments, we 
applied the following magnitude cutoffs for our "radius sam- 
ple": F435W < 22.6, F555W < 22.3 and F814W < 22.7. 
The first two limits were chosen according to Fig. [5] With 
these limits for F435W and F555W, more than 68% of the 
artificial clusters with input radii up to 9 pc that were re- 
trieved, had a measured radius better than 20%. In Fig.|4]we 
see that these limits are brighter than the 90% completeness 
limits in the high background region for sources up to ~8 pc. 
Since we do not use the radii measurements in F814W, in 
this passband the 90% completeness limit for a 3 pc source 
was applied. For the selection of the larger "resolved sample" 
we used the 90% completeness limits in all three passbands: 
F435W < 24.2, F555W < 23.8 and F814W < 22.7. 

4. We only selected sources that had an absolute difference in 
R e g between F435W and F555W of less than 2 pc. Tests re- 
vealed that mostly all sources in low, homogeneous back- 
ground regions already fulfilled this criterion and a check 
by eye showed that sources which did not pass this criterion 
were practically all contaminants due to a highly varying lo- 
cal background. We did not apply this size difference crite- 



rion for our resolved sample, since sources in this sample do 
not necessarily have accurate radii measured in both F435W 
and F555W. 

5. We used the distance of every source to its nearest neigh- 
bouring source as a criterion to filter out unreliable fits due 
to blending or crowding. A source was rejected if it had 
a neighbour within 5 pixels (the fitting radius of I shape). 
Of course this method only works when both neighbouring 
sources are in the original source list. 

6. To select out the remaining blended sources and crowded 
regions that were initially detected as a single source, we in- 
spected every remaining source in our sample by eye. We 
realize that this introduces some amount of subjectivity into 
our sample, but we are dealing with a face-on spiral galaxy 
with a high degree of irregularities in the background light 
due to the spiral structure and many crowded star forming re- 
gions. Therefore, visual inspection was unavoidable for our 
purpose of selecting a cluster sample with accurate radii. 
We created small images in all three filters for all the sources 
that fulfilled the above mentioned criteria, and by visual in- 
spection we looked for: 

- the presence of a second, separate peak within a distance 
of ~5 pixels 

- a very irregular shape 

- very small fitted aspect ratio for sources that appear fairly 
circular 

- crowded regions 

- very steep gradients in the background light that likely 
influenced the fitted radius 

- background galaxies 

When, based on one of these points, the source was not a 
clear cluster candidate with an accurately determined radius, 
we rejected the source from our sample. In Fig. [7] we show a 
number of sources that were rejected together with the main 
reason. For comparison we show a number of accepted clus- 
ter candidates in Fig. [8] Visual inspection removed another 
24, 21 and 22% of the sources from our sample in the low, 
intermediate and high background region, respectively. 

In Table [2] we summarize the results of our sample selection 
for the radius sample. Our final sample of stellar cluster candi- 
dates with accurate radii consists of 1284 sources, of which 57% 
are located in the high background region and 25% in the low 
background region. From here on we will refer to this sample of 
stellar cluster candidates simply as "clusters". The larger sample 
of resolved sources, which does not satisfy the radius difference 
criterion (criterion 4) and which magnitude cutoffs follow the 
90% completeness limits (criterium 3), consists of 7698, 6846 
and 5024 sources in F435W, F555W and F81 4W, respectively . 
This sample will be studied in a different paper (lHaas et alj2007l 
in prep.). 

5.2. Old Globular Clusters and Faint Fuzzies 

To see if there are any possible old (> 10 Gyr) GC candidates 
in our cluster sample, we applied the colour criteria B — V > 0.5 
and V-I > 0.8 to our cluster sample, typical for old MW GCs. 
There are only 70 clusters satisfying these criteria, showing that 
the majority of our cluster sample consists of young clusters, but 
a small fraction of ~5% is probably part of an old GC popula- 
tion or highly reddened. A more detailed identification of GC 
candidates in M51 will be carried out in future studies. 

We note t h at our dataset covers the same field used by 
iHwang & Leel (120061 from here on referred to as "HL06"), 
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Fig. 7. A selection of F435W images centered on sources 
that were rejected from the sample after visual inspection. 
The size of each image is 60x60 pixels or 3"x3". Rejection 
was based on the following criteria: a. neighbour too close 
by, b. likely blended, c & d. clearly blended, e. too irregu- 
lar, f. crowded, g. interfering background, h. radius fit on the 
background and i. likely a background galaxy. 
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Fig. 8. A selection of isolated stellar cluster candidates in 
F435W that passed the visual inspection. The size of each 
image is 60x60 pixels or 3"x3". 



who detect 49 "faint fuzzy" star clusters around the compan- 
ion of M5 1 , NGC 5195. Faint fuzzy clusters were introduced by 
lLarsen & Brodid J2000h . because they formed a sub-group in a 
radius-colour diagram of star clusters in NGC 1023. In Fig.|9]we 
show R e g versus V-I for the clusters in our sample. Six clus- 
ters in our sample satisfy the criteria of a faint fuzzy, namely 
0.6 < B-V < 1.1, 1.0 < V-I < 1.5 and R eS > 7 pc. These 



Table 2. The selection of the sample of stellar cluster candidates 
with accurate radii. The numbers indicate the amount of sources 
that satisfy the mentioned criterion and all the criteria mentioned 
above it, split in the low, intermediate and high background re- 
gion. 



Criterion 




Nr. remaining 






low 


inter 


high 


All sources 


35 980 


15 809 


23 647 


R cf f > 0.5 pc 


11064 


4799 


7238 


x 2 lxl < i 


10715 


4661 


7028 


F435W < 22.5 mag 








F555W < 22.3 mag 


472 


346 


1301 


F814W < 22.7 mag 








\ R F435W ~ R F555W\ < 2 P c 


457 


327 


1068 


No neighbour within 5 pixels 


418 


302 


939 


After visual inspection 


317 


239 


728 


Total sample: 1284 



25 ~ 
20 - 




-0.5 0.0 0.5 1.0 1.5 

V-I 



Fig. 9. R e ff versus V-I colour for the total sample of 1284 clus- 
ters. The upper-right corner of the dashed lines indicate the re- 
gion of the sub-group of faint fuzz y candidates (filled c ircles). 
These cuts are the same as used by iHwang & Led d2006l) to se- 
lect faint fuzzy clusters in NGC 5195. 



faint fuzzy candidates are indicated in Fig. [9] The figure shows 
that the faint fuzzy candidates indeed form a separate group in a 
radius-colour diagram and are not simply the largest clusters in 
the tail of a continuous radius distribution. 

The six faint fuzzy candidates seem randomly located in (or 
projected onto) the disc of M51. None of the 49 faint fuzzy can- 
didates of HL06 are therefore recovered in our cluster sample. 
This is because all these 49 candidates are fainter than the mag- 
nitude limits we apply. This is expected, since we apply rather 
conservative magnitude limits in order to achieve accurate radius 
measurements, also for large clusters in high background regions 
(§ 15. It . If we would drop our conservative magnitude limits from 
the selection criteria, we would have 37 of the 49 faint fuzzy can- 
didates from HL06 in our sample. The 12 remaining candidates 
are removed from our sample based on "inaccurate radii" crite- 
ria (large radius difference between F435W and F555W or ax 2 
larger than;^). The six faint fuzzy candidates in our sample are 
not in the sample of HL06, because these authors focused on the 
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Fig. 10. Cluster "212995", the largest cluster in our sample with 
a R e g = 21.6 pc (centre of the inset), and its location in M51 
(north is up, east is to the right). The projected galactocentric dis- 
tance is 1 .02 kpc. The scale of the inset image is ~ 60 x 60 pixels 
or ~ 122 x 122 pc. The scale of the big image is ~ 16 x 16 kpc. 



region around NGC 5195 and were therefore not covering the 
disc of M51. 

5.3. The largest cluster: 212995 

One cluster candidate in our radius sample, with our ID number 
"212995", clearly stands out from the other cluster candidates in 
radius. The cluster candidate, positioned at RA=13 h 29 m 51 ! i94, 
Dec=+47°11'19'.'63 and shown in Fig.[l0] has an (ellipticity cor- 
rected) R e g = 21.6 pc in F435W. The projected galactocentric 
distance of this cluster candidate is 1.02 kpc. Its brightness in 
F435W, F555W and F814W is 22.27, 21.52 and 20.48 mag, re- 
spectively, with corresponding B — V and V-I of 0.75 and 1.05, 
respectively. These colours make this source both an old GC and 
faint fuzzy candidate. Assuming the source is a cluster, we can 
make an estimate of its age and mass by using GALEV SSP 
models. However, since we only have photometry in three fil- 
ters, this estimate suffers from an age-extinction-metallicity de- 
generacy, introducing rather large uncertainties. Assuming a low 
extinction (E(B - V) < 0.1) and a metallicity in the range of 
0.2-1.0 Zq, the best estimate for the age is 2.7*H Gyr. The cor- 
responding best estimate for the mass is 2.5*^x10 Mq, which 
is a lower limit due to the underestimation of the aperture cor- 
rection for the photometry for such a large source. Assuming the 
metallicity is 0.02 Zq, the best estimates for the age and mass 
are 13. Gyr and 8.0^xl0 5 M©, respectively. However, the 
possibility of the source being a highly reddened young cluster 
is not fully excluded. There is also a possibility that this cluster 
candidate is actually a background galaxy, but this possibility is 
considered to be unlikely, since the cluster is located (in projec- 
tion) very close to the centre of M5 1 , where the extinction of the 
background source and the surface brightness of the foreground 
(M51)ishigh. 
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Fig. 11. The differences between the colours of clusters matched 
between our ACS dataset and the WFPC2 dataset of Bastian et 
al. (2005), versus the colour in the ACS dataset. The dashed lines 
show a difference of as a reference for the eye. 

6. A comparison between ACS and WFPC2 

As a test of the robustness of our methods, we compared the 
photometry and the radii of our clusters to the data of B05. 
B05 derived ages, masses, extinctions and effective radius es- 
timates of stellar clusters covering the central ~ 7.3 x 8.1 kpc 
of M51 using HST NICMOS and WFPC2 data. We calculated 
the transformation between a mosaic of the F555W WFPC2 
images and the ACS mosaic image with the GEOMAP task in 
IRAF, by identifying 10 sources by eye that were clearly visi- 
ble in both data sets. We then transformed the coordinates of all 
the clusters in their sample to our ACS mosaic frame using the 
GEOXYTRAN task in IRAF, and we matched sources which 
were within 1 pixels from a cluster in our sample, which had 
photometry in F435W (« B), F555W (~ V) and F814W (« I) in 
both data sets and which had a measured radius on the WFPC2 
data < R e ff < 12 pc. This resulted in 271 matched clusters, of 
which a few will be mismatched due to the uncertainties in the 
transformation and most importantly, geometric distortion of the 
WFPC2 images. 

For these 271 clusters we compared the B—V and V-I 
colours of B05 with our results, after we removed our filter de- 
pendent infinity corrections and Galactic foreground extinction 
corrections, since these were constant for the photometry of B05. 
In Fig. QT|we plot the difference between the ACS and WFPC2 
colours versus the ACS colour. For the mean differences we find: 

A(B-V) = (B-V) A cs - (B-V)wfpci = -0.005 + 0.005 
A(V-I) = (V-Pjacs - (V-I) WFPC 2 = 0.064 + 0.006, 

where the errors are the standard errors of the means, not to 
be confused with the standard deviations, which are 0.076 and 
0.098, respectively. This shows that our colours are in good 
agreement with B05 and that we can adopt the masses derived 
by B05 for these 271 clusters to study the mass-radius relation 
with the higher resolution of the ACS data in § 18.21 

We also compared the effective radii of the 271 matched 
clusters on the F555W image. In Fig. [12] we show the difference 
between the ACS and WFPC2 radius versus the ACS radius. No 
clear trend is visible, except that the average ACS radius of the 
clusters is slightly smaller than the WFPC2 radius. The mean 
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Fig. 12. The differences between the effective radii of the clus- 
ters matched between our ACS dataset and the WFPC2 dataset 
of Bastian et al. (2005), versus the effective radius in the ACS 
dataset. 



difference between the ACS and WFPC2 radius is 

ReffACS ~ ReH.WFPCl = ~0.34 + 0.06 pC, 

where the error is the standard error of the mean. The standard 
deviation around the mean is 0.94 pc. We note that some of 
the differences between ACS and WFPC2 radii are expected to 
be caused by contaminants in the matching procedure, as well 
as resolution effects (blends in the WFPC2 data) and a differ- 
ent treatment of ellipticity for both data sets. For the WFPC2 
data circular cluster profiles were assumed, while for the ACS 
data we used elliptical profiles with a transformation to a single 
R e s. Overall, the mean difference between the ACS and WFPC2 
radii is within the expected accuracy of the radius measurements 
(~ 0.5 pc), and Fig. [12] suggests that there are no strong radius 
dependent biases in our methods. 



7. The radius distribution 

Now that we have selected a sample of clusters with accurate 
radii, we will study the distribution of their radii and possible 
dependencies between radius, background region and colour in 
this section. Possible correlations between radius and luminos- 
ity, mass and galactocentric distance will be the subject of the 
next section. 

We show the effective radius distributiorQ of our sample of 
1284 clusters with linear bins in Fig. [T3l and with logarithmic 
bins in Fig. [14] In both figures we plot the radius distributions 
for cluster in the low and high background region separately. We 
see that the radius distribution peaks around 1-2 pc and then 
drops to a maximum radius of ~10 pc. In the remainder of this 
section we will first focus on the slope of the radius distribution 
at R e g > 3 pc and then focus on the location of the peak. 



3 Strictly speaking, the term distribution refers to linear intervals, 
i.e. N(R)dR, and the term function refers to logarithmic intervals, 
i.e. N(R)dlogR. In this work, however, we will not make this distinc- 
tion. We will only use the term distribution and we will specify the type 
of interval used when necessary. 



Fig. 13. The effective radius distribution of the total sample of 
1284 clusters, measured on the F435W image, using linear ra- 
dius bins (solid line). Also shown are the radius distribution of 
only the sources in the low background region (dotted line) and 
the sources in the high background region (dashed line). For bet- 
ter readability, only the (Poissonian) error bars of the total sam- 
ple are shown. 
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Fig. 14. The effective radius distribution of the total sample of 
1284 clusters, measured on the F435W image, using logarithmic 
radius bins (solid line). Also shown are the radius distribution 
of only the sources in the low background region (dotted line) 
and the sources in the high background region (dashed line). For 
better readability, only the (Poissonian) error bars of the total 
sample are shown. For comparison we also show a power law 
with a slope of -2.3, typical for a fractal radius distribution of 
star- forming gas clouds, and a slope of - 1 .2 as determined by 
B05 for the radius distribution of 407 star clusters in M5 1 with 
WFPC2 data. 



7.1. The slope of the radius distribution 
that the 
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power laws of the form N(M)dM oc M~ a dM, with the index a 
in the range of 1.5-2.0. Star-forming clouds also show a power- 
law radius distribution of the form N(R)dlogR oc R' 23 d log R 
dHarris & Pudritzlll994t lElmegreen & Falgarondl 19961) . For the 
clouds, the power-law mass and radius distributions are con- 
sistent with the clo uds having a fractal structure with a fractal 
dimen sion of 2.3 (Mandelbrotl ll983t [Elmegreen & Elmegreen 
2001 ). Since the mass distribution of clusters is similar to 
the mass distribution of clouds, one might naively expect the 
radius distributions also to be similar. We plotted the radius 
distribution of star clusters in M51 in Fig.[14l using logarithmic 
bins. In this figure a power law would be a straight line. We 
plotted two lines with a slope of -2.3 and -1.2. The first 
slope is consistent with the power-law distribution of the form 
N(R)dlogR oc R 23 dlogR or N(R)dR oc R~ 3 - 3 dR, observed 
for star-forming gas clouds on every length scale (down to the 
smallest scales of ~ 0.1 pc). We see in Fig. [14] that the radius 
distribution of the clusters in M51 between ~3 and ~10 pc can 
not be approximated by the same power law as the one for the 
star-forming gas cloudsQ 

The slope of —1.2, indicated in Fig. [14] is consistent with 
the power-law distribution of N(R)dR oc R~ Z2 dR, found by B05 
for 407 clusters between 2-15 pc in M51 using WFPC2 data. 
Although the slope of our observed radius distribution in the 
range ~ 3-6 pc is similar to the slope observed by B05, our ob- 
served radius distribution is considerably steeper at larger radii. 
However, we note that we used a larger sample of clusters, mea- 
sured at twice the resolution and which was checked by visual 
inspection for contaminants and blends. The cluster sample of 
B05 is therefore expected to have a larger fraction of contami- 
nants and blends than the current sample. We note however, that 
the current sample is still biased against clusters in crowded re- 
gions, but for the remainder of this study we assume that the 
current sample is representative of the entire star cluster popula- 
tion of M51. 

Fig. [14] shows that the radius distribution of star clusters in 
M51 is not consistent with a fractal structure. This suggests that 
after the formation of the clusters from a fractal gas, their radii 
have changed in a non-uniform way. Unfortunately, at the small 
radius end of the distribution a possible bias against small clus- 
ters can not be completely ruled out, since in a fractal gas the 
smallest clusters are expected to form in groups closest together. 
These small clusters could have been rejected from our sample 
by the close neighbour criterion (§ [3). Without this criterion, 
however, our sample would also be biased against small clusters 
due to blends. This bias is not expected at the large radius end of 
the distribution, where the radius distribution of star clusters in 
M5 1 is clearly not consistent with a fractal distribution. 

The picture in which the radii of the clusters change shortly 
after their formation is consistent with various suggested expla- 
nations for the lack of the mass-radius relation of clusters (see 
§ 18. 2| >. One possible explanation is that interactions between 
young star clusters and gas clouds lead t o dynamical heatin g 
and therefore expansion of the clusters (e.g. Gieles et al. 2006b). 
This expansion will be strongest for the largest and least concen- 
trated clusters because of their lower density and it will therefore 
make an initial power-law distribution more shallow. Although 
cluster-cloud interactions are probably part of the explanation 
for the lacking mass-radius relation, on its own this scenario fails 
to explain the expansion of the smallest clusters, which is neces- 
sary to weaken the mass-radius relation. 



4 Between 0.5 and ~10 pc a log-normal distribution provides a rea- 
sonable fit to the data (not shown here). 



Another suggested explanation for the weak or missing 
mass-radius relation of clusters is a star form ation efficiency 
(SFE ) which scales with the mass of the clouds dAshman & Zepfl 
200 1 ) combined with the early res i dual gas loss from cluster s 
(iHillsll 19801: iGever&BurkerthOQlt iGoodwin & Bastianll2006h . 
In this scenario the forming clusters initially follow the same 
mass-radius relation as the clouds. However, the removal of 
binding energy will lead to the expansion of preferably small 
clusters, since they form from low-mass clouds which lose rela- 
tively more residual gas. On its own, however, this explanation 
will lead to a steeper radius distribution of clusters compared to 
clouds, i.e. with a slope < -2.3, contrary to what we observe 

(Fig. ia. 

Although Fig. [14] suggests that the radii of the clusters have 
changed shortly after their formation, our radius distribution is 
consistent with neither scenario. Perhaps a scenario including a 
combination of stochastic cluster-cloud interactions, expansion 
due to residual gas expulsion and a mass-dependent SFE can 
change the radius distribution in a way that is more consistent 
with the observed radius distribution. However, a fundamental 
problem of the missing mass-radius relation of clusters is that 
there are also high-mass clusters with small radii. The existence 
of these clusters can not be explained by the scenarios mentioned 
so far, which all rely on the expansion of clusters. Therefore, 
we need better scenarios and more insight in cluster formation 
theory to account for the differences in the radius distributions 
between clouds and clusters. 

7.2. The peak of the radius distribution 

In Fig. [13] we see that there is a peak in the radius distribu- 
tion around ~1.5 pc. If we assume that the star clusters in M51 
formed from a fractal gas, this is consistent with the suggestion 
of the expansion of preferably the smallest clusters, i.e. cluster 
with R e ff < 1.5 pc which expanded to radii >1.5 pc. 

Fig. [T3land[T4l also show that the radius distribution of clus- 
ters in the low background region (the dotted lines) shows a more 
pronounced turnover, i.e. there are relatively fewer small clus- 
ters in the low background region compared to the high back- 
ground region. This suggests that generally the smallest clusters 
are mainly found in the high background regions (e.g. inside the 
spiral arms). The medians also show this trend: while the median 
R e g of our complete cluster sample is 2.1 pc, it is 1.9 pc for the 
high background and 2.7 pc for the low background region. 

We stress that it is not very likely that this trend is bi- 
ased due to selection effects, since we applied magnitude limits 
that are brighter than the 90% completeness limits in the high 
background region (§ and visual inspection removed most 
background galaxies in the low background region and possible 
blends in the high background region. Also, the more compact 
clusters are easier to detect, so it is not likely that there is a selec- 
tion effect against small clusters in the low background regions. 

Fig. [TBI shows the B-V colour distribution of clusters in the 
low and high background region. The high background region 
has a higher fraction of blue clusters (B-V < 0.05) than the low 
background region. This fraction is expected to be stronger when 
extinction is taken into account, since clusters in the high back- 
ground region are likely more strongly reddened than clusters in 
the low background region. 

Fig- 03] and Q3J suggest that there is a relation between 7? e ff 
and colour. In Fig.[16]we show the radius distribution for 2 equal- 
sized samples with B-V < 0.1 ("blue") and B-V > 0.1 ("red"). 
We indeed see a shift in the radius distribution towards larger 
radii for our red subsample. The median value follows this shift: 
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Fig. 15. The normalized B-V colour distribution for clusters in 
the low (dotted) and high (dashed) background region. The high 
background region has a higher fraction of blue clusters than the 
low background region. 



Fig. 17. The median R e g versus colour for bins containing an 
equal number of clusters for the complete cluster sample (top) 
and the low (middle) and high (bottom) background region. The 
vertical error bars indicate the range in which 68% of the clusters 
is contained. 
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Fig. 16. The effective radius distribution split in a sample with 
B-V < 0.1 (blue, dashed line, 645 clusters) and B-V > 0.1 
(red, dotted line, 639 clusters). The radius distribution of the red 
sample is shifted towards larger radii. 



for the blue sample the median R e ft is 1.8 pc, while for the red 
sample it is 2.5 pc. In Fig. [T7] we show the median 7? e ff versus 
B-V and V-I colour. Because all bins contain an equal number 
of clusters, size-of-sample effects are excluded. Both for low and 
high background regions and B-V and V-I colours we see a 
similar trend of a median R e ff increasing with colour, although 
the scatter is high and the trend is strongest for B—V colours. 

Therefore, the observed difference in the radius distribution 
between low and high background regions can be explained by a 
higher fraction of red clusters in low background regions, which 
are generally slightly larger. For young clusters, colours become 
redder with age. This is consistent with a larger fraction of blue 
clusters in high background regions, since these regions follow 
the high density spiral arms, where most clusters are expected to 
form. If the observed spread in colour is also a spread in the age 
of the clusters, the slight increase in median R e ff with colour in 



Fig. [T7] suggests a dynamical evolution of the clusters with age. 
The fact that the increase in radius is strongest for B-V colours 
supports this suggestion, because B-V is more sensitive to age 
than V-I. 

In this hypothesis, newly formed clusters in spiral arms are 
generally small, reflecting the high pressure and density of their 
parental gas clouds. In the subsequent early evolution of the clus- 
ters an increase in size is expected, likely due to dynamical heat- 
ing from cluster-cluster and cluster-cloud encounters and due 
to the removal of binding energy whe n the clusters lose mass 
dGoodwinlll997l : iBoilv & Kro upa 2003). Clusters also expand 
when moving out of th e spiral arm, due to tidal forces from 
the spiral density wave dGieles et alj |2~007). This hypothesis is 
consistent with the low background regions containing a rela- 
tively larger fraction of older, more evolved clusters with there- 
fore slightly larger radii Q 

If clusters expand, they wi ll do this on a dyna mical (cross- 
ing) timescale of a few Myrs dLada & Ladall2003h . The typical 
timescale for a cluster to move from the high to low background 
region will be about half the time between 2 spiral arm passages, 
which for a 2-armed spiral galaxy is 



t(Rc) 



ttRg 



2(v imc - £y? G ) 



(4) 



in which Vdi sc is the circular velocity in the disc and £2 p is the 
angular pattern speed. F or M5 1 this gives t( 1 kpc) « 10 Myr 
(using V = 200 km s" 1 dGarcfa-Burillo et all [1991 and Q p = 
37 kms^'kpc 1 (Zimm er et al.ll2004T) ). This is a typical lower- 
limit for the timescale to move from the high to low background 
region. For the more average galactocentric distance of 5 kpc, 



5 If one would actually quantify any age-radius relation one needs 
to be aware of possible biases, due to a slight mass-radius relation or 
size-of-sample effects. E.g. at older ages, the low-mass clusters will first 
fade below the detection limit, so any observed age-radius relation could 
then result from a possible mass-radius relation. Also, if one would let 
the absolute age intervals increase with age (i.e. logarithmic binning), 
one would sample the radius distribution up to larger radii for older ages 
and the average radius would seem to increase with age. 
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Fig. 18. The radius distribution of the 27 1 matched clusters for 
which we have age estimates from B05, split in a sample with 
log(age)>7.5 and log(age)<7.5. 

t(5 kpc) as 5.2 ■ 10 8 yr. These timescales are longer than the 
expansion timescale of the clusters, and are therefore consis- 
tent with the low background region containing a considerable 
number of larger clusters than the high background region, if the 
clusters started expanding shortly after their formation in a spiral 



7.3. A radius-age relation? 

We have used the 271 matched clusters with an age estimate 
from B05 to search for a correlation between age and R e g . These 
clusters show a weak relation between radius and age of the form 
R e ff °c age Q ' Q8±0 03 , w i th a la rge scatter. This is consistent with 
the relation iLee et aL I (120051) observe for clusters in M5 1 using 
WFPC2 data (best fit slope of 0.06 + 0.02). Fig. [18] shows the 
radius distribution for the matched clusters, split in two equal- 
sized samples with log(age)<7.5 and log(age)>7.5. The distri- 
butions are not very different, although a weak trend is visible 
since the older sample has slightly more large clusters than the 
younger sample. However, these differences are within the sta- 
tistical errors and also a K-S test did not give a significant result 
(p = 0.126). The median R e g follows a similar trend: the me- 
dian R e g is 1.8 and 2.2 pc for the younger and older population, 
respectively. 



7.4. A comparison with other results 

We compared the mean and median R e ff of our cluster sample 
to other work, but we note that these comparisons can easily be 
biased by differences in the lower limits of the radius and other 
selection criteria. The mean and median R e ff of our total sample 
are 2.5 and 2.1 pc, respectively. The mean and median R e fi of 
the 70 GC candidates in our sample are 3.7 + 0.4 and 2.4 pc, 
respectively, where 0.4 is the standard error of the mean (the 
standard deviation is 3.4 pc). If we restrict ourselves to clusters 
in the low background regions, the mean and median R e ff are 2.9 
and 2.7 pc, respectively . This last value is the same as the mean 
j?pff lJordan et alj ([2005) found for thousands of GCs obser ved in 
100 e arly-type galaxies of the ACS Virgo Cluster S«rvey. lLarsenl 
(2004) studied the effective radii of stellar clusters in 1 8 nearby 
spiral galaxies using HST WFPC2 images, and he found a mean 



R eS of 3.94 + 0. 12 pc. lLee et all 12005) studied the radii of clus- 
ters in M5 1 using HST WFPC2 observations covering parts of 
the mosaic image used in our study, and they found a mean and 
median R e g of 3.7 and 3.1 pc, respectively. The ACS camera has 
about twice the resolution of the WFPC2 camera and can there- 
fore detect small clusters more efficiently. This could explain 
why our mean and median R e s are smaller than the values from 
ILee et al] (120051) . The median R e ff for Gal actic GCs is ~3.3 pc 
dHarrislll996l) . very similar to the value iBarmbv et al] {2006) 
found for their sample of blue clusters in M101, namely 3.2 pc. 
We see that the mean and median values of R e ff of our clus- 
ter sample are somewhat smaller than what is generally found, 
although the differences become smaller when we restrict our- 
selves to the clusters in the low background regions or the GC 
candidates. 



8. Correlations between mass, radius and distance 

In this section we will combine the effective radii of the clusters 
with other parameters, such as their galactocentric distance, lu- 
minosity and for some clusters their mass. Correlations between 
these parameters for clusters in M5 1 have already been studied 
by B05, using WFPC2 data of the inner 5 kpc of the disc of M5 1 . 
We again search for correlations, but now using the ACS data out 
to a galactocentric distance of ~10 kpc. 

8. 1 . Galactocentric distance 

For Galactic GCs there is a relation between the distance to the 
Galactic centre, Rq, and the half-light diameter Z?o.5> of the form 
Dq.5 °c Rq 2 dVan den Bergh et alj[l99ll) . This same tren d is ob- 
served for the cor e radii of old clusters in the LMC (Hodge 
fl96l iMatec 



(Hesseret 



eo 
al. 



1987) and for the sizes of old GCs in NGC 5128 



1984). However, these cluster populations are all 



old and mainly residing in the outer regions or halo of their host 
galaxies, while we are studying a population of mostly young 
clusters in a disc. 

For the GCs, different explanations for the observed re- 
lation between radius and galactocentric distance have been 
suggested. One possibility could be that the sizes of GCs re- 
flect the densities of the gas clouds from which they formed, 
i.e. compact GCs preferentially formed from dense gas clouds 
near the centres of galaxies, while larger GCs preferentially 
formed in the less dense hal o regions dVan den Bergh et all 
Il991h. lHarris & Pudritzl (|1994|) use the Ebert-Bonnor relations 
(lEbertlll 955b iBonnorll 19561) to show that 



Rgmc oc m^ c p; 1/4 , 



(5) 



in which Rgmc, Mcmc and P s are the radius, mass and surface 
pressure of the gas clouds, respectively. The Van den Bergh rela- 
tion then arises naturally if the clusters form from gas clouds 
of which the surface pressure scales with the ISM pressure 
(P s « 5-1QPjsm), which in turn scales as 



ISM 



ocR 



-2 



(6) 



for the halo region (lHarris & Pudritzl 1 19941) . From Eq. [5] and [6] 
and the observation that the mean GC mass does n ot change 
with galactocentric distance (Harris & Pudritzl fl 994). the Van 
den Bergh relation follows. However, models like these assume 
that the relations with surface pressure are also valid in the cloud 
cores, where the clusters actually form, but this assumption is not 
necessarily valid. 
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Fig. 19. The median R e g versus galactocentric distance, for bins 
containing an equal number of clusters. The vertical error bars 
indicate the range in which 68% of the clusters is contained. 

Another possible explanation for the observed relation be- 
tween radius and galactocentric distance for GCs is more evo- 
lutionary in nature. It assumes that the GCs have reached tidal 
equilibrium with their host galaxy. The tidal radius of a cluster 
in an external logarithmic potential field scales as: 

r, <* M c Y 3 4 /3 , (7) 

where M c j is the cluster mass dBinnev & Tremainell 19871 chap. 
7.3). Thus, when a cluster is relaxed, in tidal equilibrium with its 
host galaxy and filling its Roche lobe, its tidal radius is expected 
to scale as r t oc /?„' . We would also expect the effective radius 

to scale as R e g oc Ra , only if the density profile of the cluster 
would be constant and if the cluster is in tidal equilibrium with 
the galaxy. For young clusters in the disc, however, the validity 
of these assumptions remains to be seen. 

In Fig.[19]we show the median R e g versus the galactocentric 
distance for the clusters in M5 1 . There seems to be a slight in- 
crease in R e ff with distance, but the scatter is large (reflected in 
the large error bars). We have tried to fit a function of the form 



^eff 

1 pc 



Rg 
1 kpc 



(8) 



to the unbinned data, in which c is a constant, and for the best 
fit we find a = 0.12 + 0.02. This relation is significantly weaker 
than the predicted a = 0.67 (for GCs in tidal equilibrium) or the 
observed relation for Galactic GCs (a - 0.5). 

The radius distribution changes for different galactocentric 
distance bins, as shown in Fig. [20] The radius distribution close 
to the centre of M5 1 (distance < 3 kpc) is shifted towards smaller 
radii compared to the radius distributions at larger distances. A 
K-S test confirmed that it is unlikely (p < 0.001) that the radius 
distribution for distance < 3 kpc and > 5.5 kpc are drawn from 
the same parent distribution. 

Since we observe a relation between R^ and color (FigfTTl), 
any increase in radius with galactocentric distance could be the 
result of an increase in colour with galactocentric distance. In 
Fig. [21] we show B-V versus galactocentric distance for the 
1284 clusters that were also shown in Fig. [19] No obvious trend 
is visible, which is expected since at all galactocentric radii 
we encounter clusters in both arm and interarm regions. The 



Fig. 20. The effective radius distributions for three different 
galactocentric distance regions in M51. The solid line shows 
clusters close to the centre of M51 (distance < 3 kpc), the 
dotted line shows clusters at intermediate distances (3< dis- 
tance <5.5 kpc) and the dashed line shows clusters at large galac- 
tocentric distances (> 5.5 kpc). The changes in the distributions 
are consistent with a slight increase in R e g with distance. 




Distance (kpc) 

Fig. 21. The median B-V colour versus galactocentric distance 
of the 1284 clusters in our sample for bins containing an equal 
number of clusters. The vertical error bars indicate the range in 
which 68% of the clusters is contained. No obvious trend is vis- 
ible in the data. 



arm regions are mostly high background regions and will there- 
fore contain mostly blue clusters, while the interarm regions 
are mostly low background regions and will therefore contain 
mostly red clusters (Fig [TBI). The observed R e g oc /j0- 12± 02 j s 
therefore not likely a result of the relation between radius and 
colour. 

B05 did not find a relation between R e g and Rq for Rq < 
5 kpc. For Rq < 13 kpc we find a weak relation, but this rela- 
tion is considerably weaker than the observed relations for old 
GCs. Therefore, the clusters we observe in the disc of M51 are 
either forming under different conditions than the GCs, or the 
observed relation for GCs emerged during their longer dynam- 
ical evolution. We consider the first explanation to be the most 
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likely one, since GCs probably formed outside a spiral disc, in 
regions where the surface pressure of their parental clouds de- 
creased with distance (Eq. [6j. For clusters in spiral arms it is 
not expected that the surface pressure scales in a similar way 
with distance. Likely, the higher pressure inside s piral arms de- 
crease s less strongly with galactocentric distance. iRix & Riekd 
(1993) find that the arm/interarm density contrast for M51 in- 
creases with galactocentric distance, consistent with this picture. 
In this case a strong radius-distance correlation for the clusters 
is not expected. 

8.2. Mass-radius relation 

One of the most peculiar properties of star clusters is the lack of a 
clear relation between their mass and radius. Star clusters are be- 
lieved to form from Giant Molecular Clouds (GMCs), fo r which 
a clea r relation between mass and radius is observed. iLarsonl 
(1981) finds that the internal velocity dispersion of GMCs, <x, 
scales with their size, Rgmc, as cr oc /&?? c . Assuming the GMCs 
are in virial equilibrium, this leads to a mass-radius relation 
of the form Rgmc k -^gmc' Also assuming virial equilibrium, 
ISolomon et all (119871) find R GMC = 1 /540 • M°^ c . These ob- 
servations are consistent with GMCs having a constant surface 
density (E = M G mc/^g MC )- 

From the Ebert-Bonnor relati ons for pressure bounded , self- 
gravitating, isothermal spheres dEberj 1 9551: iBonnorl 1 956), both 
Eq.[5]as well as: 



1/2 



(9) 



can be derived, in which P s is the surface pressure 
dHarris & Pudritzll994t lAshman & Zepfl200ll) . So the observed 
mass-radius relation and constant surface d ensity for clouds are 
expected if the surface pressure is constant (Elme greenl 19891) . 

When clusters emerge from GMCs, the mass-radius rela- 
tion appears to be erased, indicating that high-mass clusters 
have higher stellar densities than low-mass clusters. A constant 

1 /3 

stellar density would predict R e ff oc M [ , which is not ob- 
served. IZepfetalJ d 19991) find R eS oc L for young clusters 
in NGC 3256, where L is the luminosity of the cluster which 
scales directly with the mass (since their cluster sample suggests 
that colour is independent of luminosity and therefore they as- 
sume that the mass-to-light ratio is mostly independent on lu- 
minosity). For c l usters in a sample of (non-interacting) spiral 
galaxies. iLarsenl d2004) finds R e ff oc M^ w±om . The effective ra- 
dius of the old Galactic GCs also do es not seem to correlate wit h 
their luminosity and thus their mass dVan den Bergh et alj 199 lh . 
Mackev & Gilmorel d2003l) report that for a sample of 53 rich 
LMC clusters, there seems to be no strong correlation between 
their mass and core radius, either. 

In Fig. [22] we show /? e ff versus magnitude in F435W for the 
M5 1 clusters, split in the low and high background region. For 
these clusters we do not have mass estimates. However, it is ex- 
pected that the age range for the largest fraction of this cluster 
sample is not very large, because most clusters are blue and lo- 
cated in the spiral arms. Many clusters are therefore expected to 
have similar mass-to-light ratios and therefore any mass-radius 
relation should also be visible as a relation between magnitude 
and radius. Fig. l22l shows that clusters in the high background 
regions show a slight trend of radius decreasing with luminos- 
ity. Clusters in the low background regions show a less obvious 
trend, although the median /? e ff of the brightest bin is consider- 
ably larger, especially compared to the brightest bin in the high 
background region. A fit on the unbinned data points of the form 
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Fig. 22. The median R e ff versus F435W magnitude for our sam- 
ple of 3 1 7 clusters in the low background region (top) and for our 
sample of 728 clusters in the high background region. The bins 
contain an equal number of clusters and the vertical error bars 
indicate the range in which 68% of the clusters is contained. 



R e s oc L x , with L the luminosity in the F435W passband, resulted 
inx = 0.15+0.02andx = -0.1 1+0.01 for the low and high back- 
ground region, respectively. We verified that applying the size- 
dependent aperture correction of Eq. Q] instead of the constant 
aperture correction for a 3 pc source (§ I2.5K would not change 
this result considerably (x = 0.13 + 0.02 and x = -0.10 + 0.01 
for the low and high background, respectively). 

It is not likely that this observed differences in power-law 
indices is a bias due to our detection limits, since we use magni- 
tude limits (§ 15. U brighter than the 90% completeness limits for 
~8 pc sources in the high background regions. Due to the lack 
of age estimates of these clusters, there is a degeneracy between 
age and mass. Therefore it is not certain if any evolution in R e s 
with luminosity is mainly caused by age effects, mass effects, 
or a combination of both. More measurements of the ages and 
masses of clusters which also have accurate radii estimates are 
necessary to break this age/mass degeneracy (e.g. through addi- 
tional [/-band imaging). 

Using the 271 clusters that were matched with the cluster 
sample of B05 and for which we therefore have mass estimates, 
we show R e ff versus mass in Fig. [23] No apparent relation is visi- 
ble. This sample is too small to make a distinction between clus- 
ters in low and high background regions, since it mostly covers 
the inner high background regions of M5 1 . 

We conclude that we do not find evidence for any direct re- 
lation between mass and radius of the clusters, although we find 
weak relations between luminosity and radius, changing with 
background region. The suggested explanations for a lacking 
mass-radius relation were already mentioned in § 17.11 but we 
stress again that none of these scenarios are currently capable 
of explaining the observed differences in the radius distributions 
between clouds and clusters. 



9. Summary and conclusions 

We have used the HSTACS mosaic image of M5 1 to detect 7698, 
6846 and 5024 stellar clusters across the spiral disc in F435W, 
F555W and F814W respectively, based on effective radius (R e ff) 
measurements. We presented the dataset and described the meth- 
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Fig. 23. R e g versus mass for the 27 1 clusters that were matched 
with the cluster sample of B05. 

ods used to select our cluster sample, including tests with artifi- 
cial clusters to show the accuracy, limits and robustness of our 
methods. We divided the data in 3 regions with respectively a 
low, intermediate and high background, where the high back- 
ground traces the spiral arms. We selected a sample of 1284 
clusters with the most accurate radius estimates to study the ra- 
dius distribution and relations between radius, mass, luminosity, 
galactocentric distance and background region. From these stud- 
ies we conclude the following: 

1 . The effective radii of the clusters are distributed between our 
fixed lower limit of 0.5 pc and ~10 pc (Fig. [TBI. The mean 
and median R e g of our accurate radii sample are 2.5 and 
2.1 pc, respectively. This is smaller than what is generally 
found for young clusters in spiral galaxies. 

2. The radius distribution of clusters in M51 can not be fitted 
with a power law similar as the one for star-forming gas 
clouds. This suggests that shortly after the formation of the 
clusters from a fractal gas, their radii have have changed in a 
non-uniform way. 

3. 70 clusters in our sample satisfy the colour criteria for being 
old GCs. These clusters are slightly larger than the rest of the 
cluster sample (median R e g = 2.4 pc). We find 6 clusters in 
our sample satisfying the criteria for being "faint fuzzy" star 
clusters projected onto the disc of M51 (Fig. [9]). 

4. The largest cluster in our sample has R e g = 21.6 pc and 
a projected galactocentric distance of 1.02 kpc (Fig. ITOb . 
Assuming low extinction and metallicity (E(B - V) < 0.1, 
Z = 0.2-1 .0 Z©), we estimate its age to be 2.7+ Gyr and its 
mass to be 2.5+j jXlO 5 M©. Assuming extremely low metal- 
licity (0.02 Z©) results in 13.2^-8 Gyr and 8.0+^xl0 5 M©, 
for its age and mass, respectively. 

5. Comparing clusters in the low and high background regions, 
we find that the high background regions, i.e. the spiral arms, 
have a higher fraction of blue clusters, consistent with the 
idea that these regions are the preferred formation sites for 
clusters (Fig. [Tsj. 

6. We detect an increase in R e g with colour, most strongly for 
B—V. Since we detect most of the redder clusters outside the 
spiral arms, the median R e g outside the spiral arms is larger 
than inside the spiral arms: 2.7 and 1.9 pc, respectively. 
The radius distribution of clusters in the low background re- 



gion also shows a more pronounced turnover around 1.5 pc 
(Fig.Qj). We speculate that if the observed spread in colour 
is also a spread in the ages of the clusters, this observation 
suggests a dynamical expansion of the clusters with age. In 
this hypothesis, newly formed clusters in spiral arms are gen- 
erally small, their radii reflecting the high surrounding pres- 
sure of the parental gas clouds. In the subsequent early evolu- 
tion of the clusters an increase in size is expected, likely due 
to dynamical heating from cluster-cluster and cluster-cloud 
encounters and due to the removal of binding energy when 
the clusters lose mass. 

We do not observe a strong correlation between R e g and 
galactocentric distance for the clusters in the disc of M51 
out to ~13 kpc. A weak trend is visible of the form R e g oc 
^o.i2±o.02^ the scatter is large. For old GCs, mainly resid- 
ing in the outer regions or halo of other galaxies, a steeper 
relation is observed, possibly caused by the decreasing pres- 
sure of their parental gas clouds with galactocentric distance. 
The weaker relation for the clusters in M51 could be ex- 
plained by the observation that most of the clusters reside 
in the spiral arms. Since the spiral arms are expected to have 
a higher pressure and they extend out to large galactocen- 
tric distances, a strong radius-distance correlation is not ex- 
pected. 

8. We do not observe a correlation between cluster mass and 
radius for the 271 clusters of which we have mass estimates. 
We find weak relations between cluster luminosity and ra- 
dius for our sample of 1284 clusters. If fitted with a power 
law of the form R e g oc L x , we find x = 0.15 + 0.02 and 
x = -0.11 ± 0.01 for the low and high background region, 
respectively (Fig. l22l . Explanations of the lack of a strong 
mass-radius or luminosity-radius relation probably need to 
be sought in the early dynamical evolution (expansion) of the 
clusters just after their formation. Current scenarios which 
focus on the expansion of clusters due to either dynamical 
heating or the removal of binding energy due to gas expul- 
sion are not consistent with the observed differences in the 
radius distributions between clouds and clusters. 

Acknowledgements. We thank Peter Anders for useful discussions, tips and 
comments. We also thank Marcelo Mora at ESO/Garching for discussions and 
for kindly providing us with the empirical PSFs. We would like to thank Max 
Mutchler, Richard Hook and Andrew Fruchter for discussions regarding the ef- 
fects of the drizzle routine on the PSF. We thank Narae Hwang for kindly pro- 
viding us the list of faint fuzzies around NGC 5195. 



References 

Anders, P., Gieles, M., & de Grijs, R. 2006, A&A, 451, 375 
Ashman, K. M. & Zepf, S. E. 2001, AJ, 122, 1888 

Barmby, P., Kuntz, K. D., Huchra, J. P., & Brodie, J. P. 2006, AJ, 132, 883 
Bastian, N., Gieles, M., Efremov, Y. N., & Lamers, H. J. G. L. M. 2005a, A&A, 
443, 79 

Bastian, N., Gieles, M., Lamers, H. J. G. L. M., Scheepmaker, R. A., & de Grijs, 

R. 2005b, A&A, 431, 905 (B05) 
Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Bik, A., Lamers, H. J. G. L. M., Bastian, N., Panagia, N, & Romaniello, M. 
2003, A&A, 397, 473 

Binney, J. & Tremaine, S. 1987, Galactic dynamics, Princeton series in astro- 
physics (Princeton University Press) 

Boily, C. M. & Kroupa, P. 2003, MNRAS, 338, 665 

Bonnor, W. B. 1956, MNRAS, 116, 351 

Brand, J. & Wouterloot, J. G. A. 1995, A&A, 303, 851 

de Grijs, R., Bastian, N., & Lamers, H. J. G. L. M. 2003, MNRAS, 340, 197 

Dowell, J. D., Buckalew, B. A., & Tan, J. C. 2006, ArXiv Astrophysics e-prints 

Ebert, R. 1955, Zeitschrift fur Astrophysik, 37, 217 

Elmegreen, B. G. 1989, ApJ, 338, 178 

Elmegreen, B. G. & Elmegreen, D. M. 2001, AJ, 121, 1507 

Elmegreen, B. G. & Falgarone, E. 1996, ApJ, 471, 816 



R. A. Scheepmaker and M. R. Haas and M. Gieles et al.: ACS imaging of star clusters in M51 



17 



Elmegreen, D. M, Kaufman, M., Elmegreen, B. G., et al. 2001, AJ, 121, 182 

Elmegreen, D. M. & Salzer, J. J. 1999, AJ, 117, 764 

Elson, R. A. W., Fall, S. M., & Freeman, K. C. 1987, ApJ, 323, 54 

Feldmeier, J. J., Ciardullo, R., & Jacoby, G. H. 1997, ApJ, 479, 231 

Fellhauer, M. & Kroupa, R 2005, ApJ, 630, 879 

Fruchter, A. S. & Hook, R. N. 2002, PASP, 1 14, 144 

Fukui, Y., Mizuno, N., Yamaguchi, R., Mizuno, A., & Onishi, T. 2001, PASJ, 53, 
L41 

Garcfa-Burillo, S., Combes, F, & Gerin, M. 1993, A&A, 274, 148 
Geyer, M. P. & Burkert, A. 2001, MNRAS, 323, 988 

Gieles, M, Athanassoula, E., & Portegies Zwart, S. F. 2007, MNRAS, 376, 809 
Gieles, M, Larsen, S. S., Scheepmaker, R. A., et al. 2006a, A&A, 446, L9 
Gieles, M, Portegies Zwart, S. F, Baumgardt, H., et al. 2006b, MNRAS, 371, 
793 

Goodwin, S. P. 1997, MNRAS, 284, 785 

Goodwin, S. P. & Bastian, N. 2006, MNRAS, 373, 752 

Haas, M. R., Gieles, M., Scheepmaker, R. A., et al. 2007, A&A (in prep.) 

Ha§egan, M., Jordan, A., Cote, P., et al. 2005, ApJ, 627, 203 

Harris, W. E. 1996, AJ, 1 12, 1487 

Harris, W. E. & Pudritz, R. E. 1994, ApJ, 429, 177 

Hesser, J. E., Harris, H. C, Van den Bergh, S., & Harris, G. L. 1984, ApJ, 276, 
491 

Hills, J. G. 1980, ApJ, 225, 986 
Hodge, P. W. 1962, PASP, 74, 248 

Holwerda, B. W. 2005, Source Extractor for Dummies v5, STScI 
Hunter, D. A., Elmegreen, B. C, Dupuy, T. J., & Mortonson, M. 2003, AJ, 126, 
1836 

Hwang, N. & Lee, M. G. 2006, ApJ, 638, L79 

Jordan, A., Cote, P., Blakeslee, J. P., et al. 2005, ApJ, 634, 1002 

King, I. 1962, AJ, 67, 471 

Kissler-Patig, M., Jordan, A., & Bastian, N. 2006, A&A, 448, 1031 
Koekemoer, A. M., Fruchter, A. S., Hook, R. N., & Hack, W. 2002, in HST 

Calibration Workshop, ed. S. Arribas, A. Koekemoer, & B. Whitmore, STScI, 

337-340 

Lada, C. J. & Lada, E. A. 2003, ARA&A, 41, 57 
Larsen, S. S. 1999, A&A Suppl. Ser., 139, 393 
Larsen, S. S. 2004, A&A, 416, 537 
Larsen, S. S. & Brodie, J. P. 2000, AJ, 120, 2938 
Larson, R. B. 1981, MNRAS, 194, 809 

Lee, M. G., Chandar, R., & Whitmore, B. C. 2005, AJ, 130, 2128 
Mackey, A. D. & Gilmore, G. F. 2003, MNRAS, 338, 85 

Mandelbrot, B. B. 1983, The fractal geometry of nature (Freeman, San 

Francisco) 
Mateo, M. 1987, ApJ, 323, L41 
Moffat, A. F. J. 1969, A&A, 3, 455 

Mutchler, M., Beckwith, S. V. W., Bond, H. E., et al. 2005, BAAS, 37 
Pavlovsky, C, Koekemoer, A., & Mack, J. 2005, ACS Data Handbook, STScI, 

Baltimore, version 4.0 
Rix, H.-W. & Rieke, M. J. 1993, ApJ, 418, 123 
Rosolowsky, E. 2005, PASP, 117, 1403 

Sanders, D. B., Scoville, N. Z., & Solomon, P. M. 1985, ApJ, 289, 373 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Schulz, J., Fritze-v. Alvensleven, U., Moller, C. S., & Fricke, K. J. 2002, A&A, 
392, 1 

Sirianni, M., Jee, M. J., Bem'tez, N, et al. 2005, PASP, 1 17, 1049 

Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 

Van den Bergh, S., Morbey, C, & Pazder, J. 1991, ApJ, 375, 594 

Williams, J. P., Blitz, L., & Stark, A. A. 1995, ApJ, 451, 252 

Zepf, S. E., Ashman, K. M., English, J., Freeman, K. C, & Sharpies, R. M. 1999, 

AJ, 118, 752 
Zhang, Q. & Fall, S. M. 1999, ApJL, 527, L81 
Zimmer, P., Rand, R. J., & McGraw, J. T. 2004, ApJ, 607, 285 



